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Abstract 

The air-shower observatory Milagro has detected a large-scale anisotropy of un- 
known origin in the flux of TeV cosmic rays. We propose that this anisotropy is caused 
by galactic magnetic fields, in particular, that it results from the combined effects of 
the regular and the turbulent (fluctuating) magnetic fields in our vicinity. Instead of 
a diffusion equation, we integrate Boltzmann's equation to show that the turbulence 
may define a preferred direction in the cosmic-ray propagation that is orthogonal to 
the local regular magnetic field. The approximate dipole anisotropy that we obtain 
explains well Milagro's data. 



1 Introduction 



High-energy cosmic rays are of great interest in astrophysics, as they provide a complemen- 
tary picture of the sky. When they are neutral particles (photons or neutrinos), they carry 
direct information from their source [H [2J [3]. Gamma rays, in particular, have revealed 
during the past 30 years a large number of astrophysical sources (quasars, pulsars, blazars) 
in our Galaxy and beyond. In contrast, when they are charged particles (protons, electrons, 
and atomic nuclei) cosmic rays lose directionality due to interactions with the /j,G magnetic 
fields that they face along their trajectory [3]. In this case, however, they bring important 
information about the environment where they have propagated. For example, the simple 
observation that Boron is abundant in cosmic rays while rare in solar system nuclei is a very 
solid hint that cosmic rays have crossed around 10 g/cm 2 of interstellar (baryonic) matter 
before they reach the Earth. 

A very remarkable feature in the proton and nuclei fluxes is its isotropy. It is thought that 
cosmic rays of energy below 10 6 GeV are mainly produced in supernova explosions, which 
are most frequent in the galactic arms. We observe, however, that they reach us equally 
from all directions. This can only be explained if their trajectories are close to the random 
walk typical of a particle in a gas, and galactic magnetic fields seem the key ingredient in 
order to justify this picture. 

Galactic magnetic fields have been extensively reviewed in the literature [SI El [TJ El [9] . It 
is known that there is an average magnetic field of order 

B 'galactic ~ 3 fJ,G (1) 

at galactic scales. This component is the background to a second component of strength 

Brandom ~ 3 5 flG (2) 

that is regular within cells of 10-100 pc but changes randomly from cell to cell. These 
magnetic fields have frozen-in field lines and are very affected by the compressions and 
expansions of the interstellar medium produced by the passage of spiral arm waves. A 10 
TeV cosmic proton would move inside a 5 fiG field with a gyroradius of 

r s = 4~ 2x 10_3 P c > ( 3 ) 
eB 

which is much smaller than the typical region of coherence. Therefore, this proton sees the 
superposition of both components clS du TQ srular magnetic field: 

B galactic B Tan d om B re g U [ ar = B . (4) 
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Notice that the determination of the galactic field using WMAP data [TOj, [TH [12] gives 
Bgaiactic- In contrast, estimates from Faraday rotations of pulsars would be sensitive to 
the same regular B that affects the cosmic proton. According to Han et al. [T3J, [8], the 
local B should be nearly contained in the galactic plane and clockwise as seen from the 
north galactic pole {i.e., following the direction of the disk rotation), although with a small 
vertical component or tilt angle. 

At these small scales the 10 TeV proton is diffused by scattering on random fluctuations 
in the magnetic field 

5B<.B . (5) 

The interaction is of resonant character, so that the particle is predominantly scattered by 
those irregularities of the magnetic field of wave number k m l/r g . Estimates from the 
standard theory of plasma turbulence [H] indicate that 5B falls as a power law for larger 
wave numbers [13], so this component is smaller than the regular B. 

In this paper we argue that the detailed observation of the TeV cosmic-ray flux obtained 
by Milagro (T5J [16] may also provide valuable information about B and SB. In particular, 
the analysis of over 10 11 air showers has produced a map of the sky showing a large-scale 
anisotropy (a north galactic deficit) of order 10 -3 . This map, which is consistent with 
previous observations [TF], HE], remains basically unexplained. Abdo et al. [151 03] have 
discussed several possible origins: 

(i) The Compton-Getting (CG) effect [19], a dipole anisotropy that arises due to the motion 
of the Solar System around the galactic center and through the cosmic ray background. The 
anisotropy observed in Milagro's map, however, cannot be fitted by the predicted CG dipole. 
In addition, the CG anisotropy should be energy independent, which does not agree with 
the data neither. 

(ii) The heliosphere magnetic field could produce anisotropies [20l [21] that can also be ruled 
out. The Larmor radius r g sets the size of the coherence cells, and for 10 TeV protons it is 
around 2 x 10 -3 pc, significantly larger than the 5 x 10~ 4 pc (100 AU) of the heliosphere. 
Moreover, as pointed out in [151 UHL the anisotropies persist at higher energies (i.e., for 
larger distance scales), supporting the hypothesis that if magnetic fields are involved they 
are extra-heliospheric. 

Here we explore the effect of the local (regular and fluctuating) magnetic fields on the 
propagation of TeV cosmic rays reaching the Earth. Most analyses model cosmic-ray prop- 
agation with a diffusion equation [U [22], [21] , assuming certain spatial distribution of sources 
and a diffusion tensor often simplified to an isotropic scalar coefficient. This provides the flux 
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Figure 1: Angles b and //, and orthogonal vectors B, T and d for Ib = 90°. 

over an extended region around the solar neighborhood. Here we intend a different approach. 
The diffusion equation derives from Boltzmann's equation, which contains more information. 
The solution of Boltzmann's equation in the vicinity of the Earth gives the statistical dis- 
tribution function f(r,p,t), a quantity related to the intensity or surface brightness used in 
astrophysics. / provides the number of cosmic rays per unit solid angle, time and surface 
from any given direction, so it can be compared with Milagro's data pixel by pixel. 

2 Cosmic-ray distribution function 

We will treat TeV cosmic rays as a fluid that microscopically interacts only with the magnetic 
fields, and our objective is to obtain the distribution function f(f,p,t) using Boltzmann's 
equation. We will take a basic cell of radius r g and will assume that the non-turbulent 
component of the fluid is stationary and homogeneous. At these relatively small distance 
(and time) scales we can also neglect cosmic-ray sources, energy loss, or collisions with 
interstellar matter. In addition, we take the cosmic rays as protons (the dominant component 
in the flux) of E = 6 TeV (the average energy in Milagro's analysis). Finally, we will assume 
that the regular magnetic field B is on the galactic plane with a galactic longitude Ib, 
although we will show that Milagro's data favors a component othogonal to this plane (as 
found in other observations [13J). In Fig. 1 we have depicted B with Ib = 90°. 

The frequency of the (b, fi) direction in the momentum of cosmic rays reaching the Earth 
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is then proportional tcQ 

f(u) = f(b,p), (6) 

where u = p/p, b is the galactic latitude, and p is the longitude relative to the direction of 
the magnetic field B. Notice that the galactic longitude of the direction defined by u is just 
I — l B + p. 

Boltzmann's equation expresses in differential form how particles move in the six-dimensional 
phase space |23j. In our case this is just 



F ■ V u f(u) = e(uxB)-V u f(u) = . 



(7) 



Now, we separate the regular and the turbulent components both in the distribution function 
and the magnetic field: 

/ - f + Sf, 

B -> B + SB . (8) 

The components SB and Sf vary randomly from one cell to another and have a vanishing 
average value, 



SB 



[Sf] = . 



(9) 



However, there may be correlations between both fluctuating quantities. In particular, we 
will assume a non-zero value of 



e (u x SB) -V u Sf\ = e u ■ [SB x V u Sf 
= e u ■ T 



Boltzmann's equation for the regular component is then 

This equation can be also written 

u ■ ( B x V u f) + u ■ f = 



(10) 



(12) 



As u is any direction, this implies B x V u f = T, i.e., the correlation T must be orthogonal 

— * —* 

to B. Taking T in the galactic plane, 



u ■ T = T cos b sin p 



(13) 



l (E/c) 2 f{u) gives the number of particles with momentum along u per unit energy, volume and solid 
angle at£»6 TeV. 
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and expressing 

(70 COS0 O/l 

with 

u b = — sin 6 cos p u^ — sin 6 sin pu r + cos &w 2 ; = — sin + cos // w r , (15) 
Eq. (|TT|) becomes 

9/ 9/ T 

— sin a — + tan b cos u 7: — h — cos 6 sin p = . (16) 
ao op B 

This equation can be solved analytically: 

/(&>/*) = fo ^1 + sino^J + /(cos 6 cos^u) , (17) 

with /o a constant that normalizes / to the number of particles per unit volume and the 
second term any arbitrary function of the variable cos b cos \l. From the direction u we 
observe cosmic rays with p = —p u\ it is straightforward to find the relation between the 
distribution function and the flux F(b, p) of particles observed at Milagro per unit area, 
time, solid angle and energy: 

F(b,p) = ^ f(-b,p + n) . (18) 

This implies 

F(b, p) = F (1 - t sin 6) + F(cos 6 cos p) , (19) 
where t = T/(f B) and F = (E/c) 2 fo. Finally, we will expand F to second order: 

F(cos6 cos/x) ~Ficos6 cos p + F 2 (cos b cosp) 2 . (20) 

The solution in terms of the galactic longitude is obtained just by expressing p = I — l B . 
Several comments are here in order. 

(i) If F\ — F2 = 0, then the solution is a dipole anisotropy, with the minimum/maximum 
in the north/south galactic poles. This dipole is then modulated by the constants Fi^, that 
introduce an anisotropy proportional to cos b cos p (i.e., the additional anisotropy coincides 
along the directions u with equal projection on B). 

(ii) The dipole anisotropy would vanish if there were no turbulence (t — 0): B implies an 
isotropy broken by the turbulence in the orthogonal plane. In contrast, the equation does not 
say anything about the direction along B. For different boundary conditions one can find 
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Figure 2: Trihedron defined by B, cf and T, and coordinate system. B is in the galactic 
plane, whereas d and T have latitude b and 90° — bo, respectively 

solutions with a forward-backward asymmetry (implying diffusion along B) or symmetric 
solutions. In particular, F\ creates an asymmetry between the (b — 0, \i — 0) and (b = 
0,/z = 180°) directions, whereas the F% contribution is symmetric. 

(in) The dominant magnetic field B, the turbulence T, and the dipole d are always 
orthogonal to each other. For B fa Bu^ the symmetry of the galactic disc could favor a 
radial turbulence, f w Tu r , like the one that we have assumed above (see Fig. lJE However, 
one can change the latitude b Q of the dipole while keeping B on the galactic plane just by 
taking the turbulence T out of the plane. In particular, the dipole will point towards the 
arbitrary direction 6 ( see Fig. 2) if 

u-T = T(smb cos b sin \i — cos bo sin b) . (21) 

The dipole solution is in that case 

F(b, /i) = F [1 — t (sin b sin b + cos fe cos b sin //)] 

+i<icos& cos fi + F 2 (cos b cos/1) 2 , (22) 

The galactic latitude Zo of the dipole is then fixed by the orientation of B in the galactic 
plane, 

l = l B + 90° . (23) 

2 Buoyancy will mainly produce ascending turbulent cells; since Coriolis forces are negligible at these small 
time scales the compression of the (frozen-in) azimuthal field lines may result into a SB also azimuthal and 
a vertical V u Sf, which imply a radial T. 
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Figure 3: Fit of the Milagro anisotropy. 

The direction of the dipole in the basis pictured in Fig. 2 is 

Ud = cos b sin l — cos b cos l u r + sin b u z 



(24) 



3 Milagro data 

Milagro data [JS] indicate a clear dipole anisotropy, with a deficit in the north galactic 
hemisphere that peaks at 5 « 10° and AR ~ 190° (i.e., b « 72° and l Q « 293°). In Fig. 3 
we plot our fit of the data (restricted to a region in the sky), which is obtained for t = 0.003 
with Fi/Fq = 0, F 2 /Fq = 0.0003 and a magnetic field B along Is = 203°. Our simple fit, an 
approximate dipole along the direction of V u / (from b , l to — b , l — ir) provides a good 
description of Milagro's anisotropy. 



The fit implies that cosmic rays move near the Earth with a mean velocity 
v /c = -— J dfl F(u) u = -0.00059 - 0.00028 u r + 0.00157 u z 



(25) 



where N = J dfl F(u) and the basis is pictured in Fig. 2. Eq. (125]) expresses the diffusion 
velocity of the fluid (the transport flux J is proportional to Nvq), and we find that it goes 
exactly in the direction of the dipole (the term F\ would change its direction but we have 
set it to zero). 

It is important to notice that the regular magnetic field B does not need to be on the 
galactic plane (our choice above), it can rotate around the dipole axis and still give the same 
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dipole solution as far as the turbulence T is rotated as well. Doing that the only changes 
would appear in the boundary conditions (Fi and F%), but the pure dipole would provide 
the simplest solution in any case. The dipole seems to point towards 

u d = -0.35 - 0.16 u r + 0.92 u z . (26) 

Therefore, we can check if this dipole observed at Milagro and the local regular magnetic 
field B (also an observational output) are perpendicular. We will consider the values of ub 
given by Han [TBI [8] . It is found that B is basically azimuthal clockwise (a pitch angle of 
either 0° or 180° depending on the definition, which changes for different authors). However, 
the observations also indicate the presence of a non null tilt angle of 6° (a vertical component 
of order 0.3 /j,G) taking the magnetic field out of the plane. We obtain an unitary vector 

u B = 0.99 + 0.00 u r + 0.10 u z , (27) 

which implies a remarkable 

u* d . U B = -0.18 . (28) 

We think that the approximate orthogonality of these two observational vectors (we obtain 
an angle of 100°) provides support to the model presented here. 

Notice that our framework could also accommodate other anisotropics in the flux, added 
to the dipole one, as far as they have the same value in all the points with equal projection 
(cos b cos /i) on B. To explain a pointlike anisotropy like the one named as region A in 
[T6] , the anisotropy itself should be along the direction of the dominant magnetic field B 
(orthogonal to d). Region A, however, is at (pA ~ —30°, I a ~ 215°), forming an angle of 58° 
with the dipole. 



4 Summary and discussion 

Although charged cosmic rays do not reveal their source, the study of their flux from different 
directions is of interest in astrophysics because it brings valuable information about the 
interstellar medium. In particular, the per mille deficit observed by Milagro could be caused 
by the local (at distances of order r g ) magnetic fields. 

Using Boltzmann's equation we have shown that the interplay between the regular and 
the turbulent components in these magnetic fields always produces a dipole anisotropy in 
the cosmic-ray flux. We find that (i) the direction of this anisotropy is orthogonal to the 
regular B and (ii) its intensity is proportional to the fluctuations 5B/B at the wave number 
k = l/r g . These two simple results have already non-trivial consequences. In particular, 
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(i) implies that a north-south galactic anisotropy would only be consistent with a dominant 
B laying in the galactic plane, whereas (ii) explains that the anisotropy is larger for more 
energetic cosmic rays: their gyroradius r g is larger, the resonant wave number k smaller, so 
the expected value of SB/B will be larger. 

We have argued that Milagro's data can be interpreted as a dipole anisotropy pointing to 
a well defined direction in the north galactic hemisphere, namely, (bo m 67°, Iq m 284°). Our 
model provides a remarkable fit of the data, so we conclude that it explains satisfactorily 
the large-scale anisotropy found by Milagro. The model implies that the dominant magnetic 
field near our position must be in the plane orthogonal to the dipole (B, the turbulence 
correlation T and d define a trihedron). This plane forms an angle 9 = 23° with the galactic 
disc. 

The data obtained by Milagro (energy, direction and nature of over 10 11 primaries) shows 
that the 10~ 3 deficit in the cosmic-ray flux from the north galactic hemisphere already seen 
in previous experiments [T71 [T5] is actually very close to a dipole anisotropy. We think 
that the analysis of the flux after substracting this dipole anisotropy could reveal further 
correlations. 
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